Optical properties of graphene 
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Reflectance and transmittance of graphene in the optical region are analyzed as a function of 
frequency, temperature, and carrier density. We show that the optical graphene properties are 
determined by the direct interband electron transitions. The real part of the dynamic conductivity 
in doped graphene at low temperatures takes the universal constant value, whereas the imaginary 
part is logarithmically divergent at the threshold of interband transitions. 

PACS numbers: 



I. INTRODUCTION 

Monolayer and bilayer graphenes are gapless two- 
dimensional (2D) semiconductors whereas its 3D prede- 
cessor, graphite, is a semimetal. Hence the dimensional- 
ity effects can be studied for the unique substance. Mono- 
layer graphene has a very simple electron band struc- 
ture. Near the energy e = 0, the energy bands are cones 
£■1,2 (p) = at the K points in the 2D Brillouin zone 
with the constant velocity parameter v = 10 8 cm/s. Such 
a degeneration is conditioned by symmetry because the 
small group C^ v of the K points has a two-dimensional 
representation. 

While the carrier concentration is decreasing in the 
field gate experiment, the graphene conductivity at low 
temperatures goes to the finite minimal values. Much 
theoretical efforts have been devoted to evaluate the 
minimal conductivity in different approaches. Theoret- 
ical and experimental researches show that the main 
mechanism of the carrier relaxation is provided by the 
charged impurities and gives the collision rate t -1 ~ 
2n 2 e i ni mp /fie 2 g £ 1 where e g is the dielectric constant of 
graphene, e is the characteristic electron energy (of the 
order of the Fermi energy or temperature), and rii mp is 
the density of charged impurities per the unit surface. 
The general quantum expression for the conductivity de- 
pending on the frequency to and wave vector k has been 
derived in our paper—. In various limiting cases, our re- 
sult coincides with the formulas of Refs^^. This expres- 
sion is valid under a restriction that the collision rate of 
carriers is less than the frequency and spatial dispersion, 
t _1 <C uj,kv. In the optical range, we can neglect the 
spatial dispersion of conductivity compared with the fre- 
quency dependence and use the dynamical conductivity 
cr(w) to study the graphene optical properties^. The op- 
tical visibility of both monolayer and bilayer graphene 
was theoretically studied in Ref. 5 focusing on the role of 
the underlying substrate. Recently, graphene transmit- 
tance spectra were observed* and the dynamical conduc- 
tivity was found as <j{uj) = e 2 /Ah in agreement with the 
theoryi 2 ^. 

In the present paper, we analyze in detail the optical 
properties of the graphene monolayer and multilayers. 
We use the dynamic conductivity a as a function of fre- 



quency u>, temperature T, and chemical potential u. The 
chemical potential of ideal pure graphene at any temper- 
ature is at the crossing of the bands E\^. With the help of 
the gate voltage, one can control the density and type (n 
or p) of carriers varying their chemical potential. Using 
the dynamic conductivity and the appropriate boundary 
conditions at interfaces, we calculate the reflection and 
transmission coefficients of graphene monolayer and mul- 
tilayers with varying carrier concentration. 



II. OPTICAL CONDUCTIVITY OF GRAPHENE 

The general expression for the conductivity used here 
is obtained in the previous paper—. For high frequencies, 
lo ^> (fcw,r _1 ), the dynamical conductivity [see Eq. (8) 
in Ref. 1 ] is given by 
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where fo{e) — {exp[(e — fj)/T] + 1} 1 is the Fermi func- 
tion. 

The first term in Eq. fl} corresponds to the intra- 
band electron-photon scattering processes. Integrating, 
we obtain explicitly: 
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^m[2cosh(/V2T)]. (2) 



Here, we write u> + ir^ 1 instead of u> in order to take into 
account the electron - impurity scattering processes. In 
such a form the intraband conductivity coincides with 
the Boltzmann-Drude expression. For the Fermi-Dirac 
statistics, u 3> T, the intraband conductivity takes the 
form 



ie 2 \n\ 



Trh(uj + it 1 ) 



(3) 



where the chemical potential determines the carrier con- 
centration no = (fi/hv) 2 /tt. 

The second term in Eq. |T]), where J — ^ is infinites- 
imal quantity determining the bypass around the inte- 
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grand pole, owes its origin to the direct interband elec- 
tron transitions. The integral is easily evaluated in the 
case of zero temperature: 



» 



e 
4h 



2 r 



V ' ' 2tt (w - 2/i) 2 



(4) 



Here and below, we consider for simplicity the case of 
the positive /i. The step function 0(u> — 2/i) conveys the 
condition for the interband electron absorption at low 
temperatures. The expression shows that the interband 
contribution plays the leading role around the absorption 
threshold w r* 2/i comparatively with the Boltzmann- 
Drude intraband term, Eqs. © and ([3]), which is impor- 
tant at relatively low frequencies, u> < fi. The logarithmic 
singularity is cut off with temperature. At the finite, but 
low temperatures, the following substitutions should be 
made in Eq. (0J 



1 



(w - 2/i) -s- - + - arctan[(w - 2/x)/2T] (5) 

2 7T 

(u) - 2/i) 2 ->• (u; - 2/i) 2 + (2T) 2 . 



The electron relaxation processes produce similar smear- 
ing. 

It is useful for numerical calculations to present the 
difference of the Fermi functions in the second integrand 
© as 



G(e) 



sinh(e/T) 



cosh(/i/r) + cosh(e/T) 



Adding and subtracting G(u>/2) in the numerator of the 
integrand, and noticing that the the principal value of 
the integral with G(u>/2) equals to zero, we arrive at the 
integral without singularities. Then we can write the in- 
terband conductivity in the form available for numerical 
calculations: 



a mter (uj) 
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(6) 

For the Fermi-Dirac and Boltzmann carrier statistics, 
correspondingly, the first term is given asymptotically 
by 



G(w/2) = 



0(w-2ii), M>T, 
tanh(w/4T), fi<T. 



(7) 



The conductivity calculated with the help of Eqs. (O, 
© is shown in Fig. [T] The step function behavior of 
the real part (absorption) and logarithmic singularity of 
the imaginary part are clearly seen at low temperatures. 
While increasing the temperature, the transition from 
the Boltzmann statistics to the Fermi-Dirac statistics is 
also evident. 

By using the gate voltage, one can control the den- 
sity of electrons (no) or holes (—no). Then the chemical 
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FIG. 1: (color online). Imaginary (left) and real (right) parts 
of conductivity in units of e 2 /h for graphene doped with the 
carrier density no = 10 10 cm -2 at temperatures noted at lines. 
The chemical potential equals to 135, 65, and 33 K at 3, 100, 
and 200 K, correspondingly. 



potential is determined by the condition 



n 



ir(hv) 2 



e[f (e-fx)-f (e + Li)}de, (8) 
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From this expression and Fig. [2ja), one can see that the 
chemical potential goes to zero while the temperature 
increases. 
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FIG. 2: (color online), (a) Chemical potential (in K) as a 
function of temperature at carrier densities noted at curves, 
(b) Multilayers sample and geometry of wave scattering. 
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III. SPECTROSCOPY OF GRAPHENE LAYERS 

In order to calculate the graphene reflectance, we apply 
Maxwell's equations 



V(V • E) - V^E 



c z 



(9) 



where eo is the ion contribution into the dielectric con- 
stant and j is the conductivity current. We consider the 
case of the p-polarization, when the field E lies in the xz 
plane and the current j has only the in-layer x component 
(see Fig 12b). 



A. Optics of a monolayer. 
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Consider the graphene monolayer at z = with eo = e g 
deposited on the substrate (z > 0) with the dielectric 
constant eo = e s (for graphene suspended in the vacuum, 
e s = 1). The ac field is given by the sum of incident and 
reflected waves in the vacuum, z < 0, and by the trans- 
mitted wave in the substrate. In the geometry consid- 
ered, the current in graphene monolayer can be written 
in the form 



j x = a(uj)S(z)E x 



(10) 



Making use of the Fourier transformations with respect 
to the x coordinate, E oc e lkxX , we rewrite the Maxwell 
equations © as follows 



dE z d E x 

dz dz 2 
dE x , 2 
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c 2 
c z 



-3a 



(11) 



0. 



The boundary conditions for these equations at z = 
are the continuity of the field component E x and the 
jump of the electric-induction z component at the sides 
of the monolayer: 



FIG. 3: (color online). Transmittance spectrum of graphene 
with carrier densities no = 10 10 cm -2 (left) and no = 10 11 
cm -2 (right) versus the frequency at temperatures noted at 
the curves; normal incidence. For the carrier density no = 
10 11 cm -2 , the chemical potential equals to 428, 389, and 294 
K at 3, 100, and 200 K, correspondingly. 



Using the boundary conditions, we find the reflection 
(r) and transmission (t) amplitudes 

1 " C ■ 2 (M) 



1 + C 



t = 



1 + C" 



where C — kl[(4ira(ui)/uj) + (e s /k s )] . The very simple 
results can be written for graphene suspended in the vac- 
uum, when e s = 1 and k s = k\. Then, the coefficient C 
is close to unity, 



C = 1 + 4ttct(uj) cos 9/c , 



(15) 



where 8 is the incidence angle. Therefore, Eqs. (fT"4|) yield 
for the reflected and transmitted amplitudes 

r = —2tt(j(uj) cos 9/c , t = 1 — 2tt(j(u) cos 9/c, (16) 



e.,-E z L=+n - £U z= _o = 4.7T 



+o 



p(ui,k x ,z)dz . (12) 



The carrier density is connected to the current in Eq. 
(|10[) according to the continuity equation 

p(u>, k x , z) = j x (oJ, k Xl z)k x /uj. 

Substituting E z from the second Eq. (fTTj) into (fT2j) . we 
find the second boundary condition 



e. s dE x 



1 dE x 



47rcr(w) 



k 2 dz 
where 

k s = yje s (u/c) 2 



dz 



E 



iu> 



The transmission coefficient \t\ 2 calculated with the 
help of Eqs. (HI),©, and (fT4"|) for normal incidence is 
shown in Figs. [3] as a function of frequency, tempera- 
ture, and carrier density. One can see that the devia- 
tion of the graphene transmittance from unity is propor- 
tional to the dimcnsionless parameter e 2 max(T, jj,)/hcuj, 
i.e., it is relatively large at the very low frequencies 
uj « e 2 max(T, p) /he, when the intraband conductivity 
plays the leading role. At higher frequencies, the inter- 
band transitions play the leading role. Then, the trans- 
mittance does not depend on frequency, being controlled 
(13) by the fine structure constant e 2 /Tic. For instance, at 
normal incidence and T <C [i < <jj/2, the transmission 
coefficient is given by 



V("/c) 2 ~ k 2 



47T 



Re <r(u>) = 1 — 7r 



he 



(17) 
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where the logarithmic terms is omitted since it con- 
tains the fine structure constant squared. As was found 
recently^, the linear (in the fine structure constant) effect 
can be observed in graphene at the visible frequencies. 

Reflectance increases with the temperature, because 
the carrier density increases. As the chemical potential 
decreases with the temperature, the temperature depen- 
dence of reflectance is not monotonic as seen clearly in 

Fig. m 



B. Spectroscopy of graphene multilayers. 

Let the multilayers cross the z axis at points z n — nd, 
where d is the distance between the layers (see Fig. [5b). 
Such a system can be considered as a model of graphite 
since the distance d = 3.35A in graphite is larger than 
the interatomic distance in the layer. So, we describe the 
carrier interaction in the presence of ac electric field with 
the help of self-consistent Maxwell's equations ©. For 
the x component of the field E x , they give 

d 2 \ 
— + k 2 s +2k s V^25{z-nd)\E x = 0, (18) 

ii ) 

where T> = 2mu(uj)k s j t g Lo . 

For the infinite number of layers in the stack, the so- 
lutions of Eq. (|18p represent two Bloch states 

e 12 ( z ) = e ±ik * nd {sin k s (z - nd) - e Tik * d x 
sin/c s [z - (n + l)d]}, nd < z < (n + l)d (19) 

with the quasi-momentum k z determined from the dis- 
persion equation 

cos k z d = cos k s d — T> sin k s d . (20) 

The dispersion equation describes the electric field ex- 
citations in the system. The quasi-momentum k z can 
be restricted to the Brillouin half-zone < k z < ir/d, 
if the parameter V is real. In the general case, while 
taking the interband absorbtion into account, we fix the 
choice of the eigen-functions in Eq. (TP9"]) by the condition 
Im k z > so that the solution e\ decreases in the positive 
direction z. 

In the long- wave limit, k Zl k s <C l/d, the dispersion 
equation (|20|) can be simplified. In this case, the di- 
electric permittivity can be introduced not only in the 
normal z direction, e zz — e g ~ 2.5, but in the tangential 
direction as well, namely, 

e xx = e g + 4:nia(uj)/ujd , (21) 

where cr{uS) is the dynamic conductivity, Eq. (fTJ), of the 
one layer. Then the dispersion equation (f2"0"l) takes the 
form 

kle xx + k z e zz = (uj/c) 2 e xx e zz . (22) 




FIG. 4: (color online). Transmittance and reflectance at the 
normal incidence for a plate with the graphene multilayers 
and the carrier density no = 10 10 cm -2 in a layer; the distance 
between layers d = 3.35A, thickness of the plate is I = lOOd, 
temperatures are noted at curves. 



We see from this equation, that the weakly damped so- 
lutions (for normal incidence, k x — 0) exist, if the real 
part of the dielectric function e xx is positive and larger 
than the imaginary part. According Eqs. and Q, 
this condition is fulfilled just below the threshold of in- 
terband transitions. Due to these weakly damped waves, 
the maximum in transmittance and the corresponding 
minimum in reflectance of the multilayer system appear 
at temperatures below the threshold as shown in Fig. [4] 

Calculations done for a thin plate with the stack of 
graphene multilayers give the following result for the 
transmission amplitude 



(fc. + A*) 2 /-(fe.-*4) a /- 1 ' K 1 

where / = exp(— ik z l), k\ = (cu/c) cos9, k z is determined 
by the dispersion equation (f2"2")l at fixed values of u> and 
k x . 

At higher frequencies, Li > 2u = 270 K (in Fig. 
SJ), after the downfall, the reflection coefficient corre- 
sponds with the interband absorption presented by the 
0— singularity in the real part of conductivity. In contrast 
to the monolayer, the effect of carriers in the multilayers 
is controlled by the large dimensionless factor (e 2 /TitgLod) 
[see Eq.f l2~Tj)] of the function varying rapidly at u> w 2u. 
It means that the interband absorption can be observed 
at relatively large frequencies and low temperatures. No- 
tice that observations of the absorbtion threshold pro- 
vide a direct method of carrier density characterization 
in graphene. 
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IV. CONCLUSION 

In conclusion, we have developed the detailed micro- 
scopic theory of the graphene monolayer and multilayers 
spectroscopy. We have shown that the reflectance from 
the monolayer is determined for infra-red region by the 
intraband Drude-Boltzmann conductivity and for higher 
frequencies by the interband absorption. We have argued 
that at low temperatures and high carrier densities, the 



reflectance from multilayers has the sharp downfall with 
the subsequent plateau. These features are caused by ex- 
citations of propagating waves and the direct interband 
electron transitions. 
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